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ABSTRACT 

This paper mainly proposes an efficient modified 
particle swarm optimization (MPSO) method, to 
identify a slider-crank mechanism driven by a field- 
oriented PM synchronous motor, dynamic 
formulations of a slider-crank mechanism have been 
successfully formulated with only one independent 
variable. The parameters of many industrial machines 
are difficult to obtain if these machines cannot be 
taken apart. In system identification, we adopt the 
MPSO method to find parameters of the slider-crank 
mechanism. This new algorithm is added with 
“distance” term in the traditional PSO’s fitness 
function to avoid converging to a local optimum. It is 
found that the MPSO method can obtain optimal high- 
quality solution, high calculation efficiency, and its 
feasibility and effectiveness. Finally, the comparisons 
of numerical simulations and experimental results 
prove that the MPSO identification method for the 
slider-crank mechanism is feasible. 

Keyword: Slider-crank mechanism, Modified particle 
searm optimization. 

I. INTRODUCTION 

A slider-crank mechanism is widely used in gasoline 
and diesel engines, and has been studied extensively 
in the past three decades. The responses of the system 
found by Viscomi and Ayre [1] are to be dependent 
upon the five parameters as the length, mass, 
damping, external piston force and frequency. The 
steady-state responses of the flexible connecting rod 
of a slider-crank mechanism with time-dependent 
boundary effect were obtained by Fung [2], A slider- 
crank mechanism with constantly rotating speed was 
controlled by Fung et al. [3]. The mathematical model 
of the coupled mechanism of a slider-crank 
mechanism was obtained by Lin et al. [4], where the 
system is actuated by a field-oriented control 
permanent magnet (PM) synchronous servomotor. 


However, the dynamic formulations of a slider-crank 
mechanism with one degree-of -freedom have more 
than one independent variable in the past researches 
[3-4], In this study, the dynamic formulation is 
expressed by only one independent variable of 
rotation angle. Moreover, its dynamic responses are 
compared well with the experimental results. 

Particle swarm optimization (PSO) is a stochastic 
population based optimization approach, and was first 
published by Kennedy and Eberhart in 1995 [5-6], 
PSO has been shown to be an efficient, roust and 
simple optimization algorithm. Most studies of the 
PSO are empirical with only a few theoretical 
analyses that concentrate on understanding particle 
trajectories. A large amount of research has been 
proposed to improve the performance of PSO. From 
these studies, much effort has been invested to obtain 
a better understanding of its convergence properties. 
These studies concentrated mostly on a better 
understanding of the basic PSO control parameters, 
namely the acceleration coefficients, inertia weight, 
velocity clamping, and swarm size [7-9]. In this 
paper, a modified PSO algorithm is proposed to 
improve the searching ability and prevent from being 
trapped in a local optimum. The main difference of 
the MPSO from the PSO is its fitness function 
considers the “distance” to avoid converging to a local 
optimum. From these empirical studies it can be 
concluded that the MPSO is sensitive to control 
parameter choices, specifically the inertia weight, 
acceleration coefficients and velocity clamping. 
However, wrong initialization of these parameters 
may lead to divergent or cyclic behavior. 

This study successfully demonstrates that the dynamic 
formulation can give a wonderful interpretation of a 
slider-crank mechanism in comparison with the 
experimental results. Furthermore, a new identified 
method using the MPSO is proposed, and it is 
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confirmed that the method can perfectly searches the 
parameters of the slider-crank mechanism driven by a 
servomotor through the numerical simulations and 
experimental results. 


II. DYNAMICS FORMULATION OF A 
SLIDER-CRANK MECHANISM 

A slider-crank mechanism is a single looped 
mechanism with a very simple construction shown in 
Fig. 1(a) and the experimental equipment of a slider- 
crank mechanism is shown in Fig. 1(b). It consists of 
three parts: a rigid disk, which is driven by a 
servomotor, a connecting rod and a slider. 

Figure 1 (a) shows the physical model of a slider-crank 
mechanism, where the mass center and the radius of 
the rigid diskare denoted as point “0” and length “ r 
”, respectively. And “/” is denoted as the length of 
the connected rod AB . The angle 6 is between OA 
and the X-axis, while the angle <j) is between the rod 
AB and the X-axis. In OXY plane, the geometric 
positions of gravity centers of rigid disk, connected 
rod, and slider, respectively, are as follows: 

Acg = Ticg=° (!) 

x 2 c g = r cos 0 + cos ^ yicg = ^ l sin (j) (2) 

A cg =rcos0+lcos</>, y 3cg = 0. (3) 

The mechanism has a constrained condition as 
follows 

r sin 6 = / sin <j>. (4) 

The angle (j) can be found from Eq. (4) as 

(j> = sin -1 sin 6 j. ) 

In the kinematic analysis, taking the first and second 
derivates of the displacement of slider B with respect 
to time, the speed and acceleration of slider B are as 
follows 

in 6 - l<fk in (j) (6) 

= -r^in O-rffecosO-lffis in^-/^cos^. (?) 


Similarly, the angular velocity ^ and acceleration^ 
are obtained as follows 




r$kosO 

lcos(f> 

rffltos <j> cos 9+r$$kos 0 sin <j>- sin 9 cos ^ 


l COS (f) 


( 8 ) 

(9) 


A machine model of a PM synchronous motor can be 
described in a rotor rotating by [10] as follows: 


T, = R J q + H + w Ai 

(10) 

1 

+ 

II 

(11) 

Where 

. ^ 

<1 

II 

(12) 

> 

f —i 

"a 

S 

+ 

. 

II 

(13) 


In the above equations, Vj and Vq are the d and q axis 
stator voltages, i d and / are the d and q axis stator 
currents, L d and L q are the d and q axis inductances, 
X d and x are the d and q axis stator flux linkages, 
and R and w s are the stator resistance and inverter 
frequency, respectively. In Eq. (13), / is the 
equivalent d -axis magnetizing current, and L md is the 
d -axis mutual inductance. The electric torque is 
*1 = | P [L md I fd i q +(L d -L q )i d i q \ (14) 

and the equation for the motor dynamics is 

7, = A + B Jd + J ,A- ( 15 ) 

In Eq (14), p is the number of pole pairs, z m is the 
load torque, B m is the damping coefficient, m r is the 
rotor speed and j m is the moment of inertia. The basic 

principle in controlling a PM synchronous motor 
drive is based on field orientation. The flux position in 
the d-q coordinates can be determined by the shaft- 
position sensor because the magnetic flux generated 
from the rotor permanent magnetic is fixed in relation 
to the rotor shaft position. In Eqs. (13-14), if i d =o the 

d -axis flux linkage x d is fixed since L wd and / are 

constant for a surface-mounted PM synchronous 
motor, and the electromagnetic torque r e is then 

proportional to ; , which is determined by closed-loop 

control. The rotor flux is produced in the d - axis 
only, and the current vector is generated in the q - axis 
for the field-oriented control. As the generated motor 
torque is linearly proportional to the q - axis current as 
the d- axis rotor flux is constant in Eq. (14), the 
maximum torque per ampere can be achieved. 

With the implementation of field-oriented control, the 
PM synchronous motor drive system can be 
simplified to a control system block diagram, as 
shown in Fig. 2, in which 
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K, =~PL md l fd , 


HM = 


l 


Js + B_ ’ 


(17) 

(18) 


where r m is the torque applied in the direction of w r , 
and the variables a> r and are the angular speed and 
acceleration of the disk, respectively. 


where /* is the torque current command. By 

substituting (16) into (15), the following applied 
torque can be obtained as follows: 

r m =K t i q -J m (Sc-B m(0r (19) 


According to Ref. [3], the dynamic equation of the 
slider-crank mechanism can be concluded to be 


M(eyifrrN{d,$)= F(e) 


( 20 ) 



(2 m 3 + m 2 )-\ —-rcos0 
c 


f 3 

r 


\ 1 


cos#sin 2 6 


+ {m 2 + m 3 )r 2 sin 2 6 + — mi — j (rcos #) 2 +— m x r 2 + J m 
3 \c) 2 


N = <m 2 r sin(9cos<9 


/ 2 r (j r f r 3 

1-- + -cos# + -^-cos 2 0 + —-cos <9 sin 2 6 

3c 2 c 3c 4 2c 3 


r 4 cos 2 #sin 2 # r 3 . 2 ^ ^ 

+-t -+ —sin <9 cos 6 


c c 

F = K t i q -(f b + F E )rsm6\ 1 + —cos 6 |. 
c - a// 2 -r 2 sin 2 6 . 


- m , — sin 3 6 + m,r sin 0 cos 6 

2 2c 3 


y 1 2^ 

1 —^-sin 2 9 + —rcos 2 9 +—cos 9 
c c c 


- m 3 —sin 3 9 + BJ§^v —m 2 gr cos 9, 
c 2 


The system becomes an initial value problem and can 
be directly integrated by using the fourth order 
Runge-Kutta method. 

III. IDENTIFICATIONS OF THE SYSTEM 
A. Particle swarm optimization 

Birds (particles) flocking optimizes a certain objective 
function in a PSO system. Each agent knows its best 
value so far (pbest) and its position. This information 
is analogy of personal experiences of each agent. 
Moreover, each agent knows the best value so far in 
the group (gbest) among pbests. This information is 
analogy of knowledge of how the other agents around 
them have performed. The PSO concept [5-6] consists 
of changing the velocity of each particle toward its 
pbest and gbest locations. In the PSO, each particle 
moves to a new position according to new velocity 
and the previous positions of the particle. This is 
compared with the best position generated by 
previous particles in the fitness function, and the best 
one is kept; so each particle accelerates in the 
direction of not only the local best solution but also 
the global best position. If a particle discovers a new 
probable solution, other particles will move closer to 
it to explore the region more completely in the 
process. 


In general, there are three attributes, current position 
Xj , current velocity v. and past best position pbest ., 

for particles in the search space to present their 
features. Each particle in the swarm is iteratively 
updated according to the aforementioned attributes. 
For example [5-9], the jth particle is represented as xj 
= (xj,l, xj,2,..., xj,g) in the g-dimensional space. The 
best previous position of the jth particle is recorded 
and represented as pbestj = (pbestj,l, pbestj,2,..., 
pbestj,g). The index of best particle among all 
particles in the group is represented by the gbestg. 
The rate of the position change (velocity) for particle j 
is represented as vj = (vj,l, vj,2,..., vj,g). The 
modified velocity and position of each particle can be 
calculated using the current velocity and distance 
from pbestj,g to gbestj,g as shown in the following 
formulas [9]: 

+ Cj RandQ • (pbestj - x (t) ) 


v ( ' +1 >=W () 

j,g j,g 


+ c 2 Rand’ () • {gbest - x\ l> ), 


, 0 ) 


J,g 


J,g j 


( 21 ) 


v (f+i) _ v (0 , V (<+1) 
x j,g ~ x j,g + v j,g ■ 

j = 1,2,A ,n;g = 1,2,A ,m 

where n is the number of particles in a group; m is the 
number of members in a particle; t is the pointer of 
iterations (generations); is the velocity of the 

particle j at iteration t, E™ 1 " < i/f < V " VdX ; w is the 
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inertia weighting factor; c,, c 2 are the acceleration 
constants; Rand(), Rand*() are random numbers 
between 0 and 1 ; v'7 is the current position of 

particle j atiteration t; pbestj is the pbest of particle j; 
gbestg is the gbest of the group g. 

In the above procedures, the parameter F" mx 

determine the resolution or fitness, with which regions 
are searched between the present position and the 

target position. If F g max is too high, particles might fly 
past good solutions. If V" mx is too low, particles may 
not explore sufficiently beyond local solutions. 

The constants c x and c 2 represent the weighting of 
the stochastic acceleration terms that pull each 
particle toward pbest and gbest positions. Low values 
allow particles to roam far from the target regions 
before being tugged back. On the other hand, high 
values result in abrupt movement toward or past target 
regions. 

Suitable selection of inertia weighting factor w 
provides a balance between global and local 
explorations, thus requiring less iteration on average 
to find a sufficiently optimal solution. As originally 
developed, w often decreases linearly from about 0.9 
to 0.4 during a run. In general, the inertia weighting 
factor w is set according to the following equation [ 6 - 

7 ]: 

w = w max - Wmax ~ Wmin X iter (22) 

*' fe r max 

where iter mdx is the maximum number of iterations 

(generations), and iter is the current number of 
iterations. 

B. Modified particle swarm optimization 

The main point of the MPSO differs from the PSO is 
to consider the “distance” in its fitness function to 
avoid converging to a local optimum. Assign a rank 
(i.e., the number place 1, 2, 3, ..., etc.)RE).tothe 

calculated error of each new individual, v k , 
k- 1,A ,PS , PS is the population size. A combined 
population with 2 x PS individuals is formed. Unlike 
previously developed statistic methods, the concept of 
“distance” is added to the fitness function to prevent 
from being trapped in a local minimum. The fitness 
score of the kth individual is modified by [ 11 - 12 ] 

F k =RE k + pxRD k ,k = l,A ,2xPS. (23) 


where p is an adaptive decay scale, p imx is set as 0.7 
and p min is set as0.005 in this paper. RD k is the rank of 
D k assigned to the kth individual, where D k is the 

distance from the individual to the current best 
solution vector, and is given by 

A=lh _v J ( 24 ) 

where v k is the vector of the kth individual in the 
combined population, and v best is the current best 
solution vector. 

An adaptive scheme is defined as [11] 

PA = R X (/Uax - Pmin ) 1 Smax ( 25 ) 

( +n = [rtg)-PA, F min (g) = F min (g- 1 ) 

W) 3 F min (g)<F min (g-l) 

(26) 

and 

p(g + 1) = p mia ; if p(g) -p A < p min (27) 

where p A is the step size; F mm is the minimum value 
of fitness functions; R is the regulating scale and is 
set as 1.25 in this paper, and g max is the maximum 
allowable number of iterations. 

Individuals will be ranked in ascending according to 
their fitness scores by a sorting algorithm. The PS 
individuals are transcribed along with their fitness for 
the next generation. If the new population does not 
include the current best solution, the best solution 
must be replaced with the last individual in the new 
population. In addition, a gradually decreased decay 
scale can satisfy a successive statistic searching 
process by first using the diversification (bigger p) to 
explore more regions, and then the intensification 
(smaller/?) to exploit the neighborhood of an elite 
solution. The current best solution (point A) for a 
minimum fitness problem as shown in the Fig. 3 may 
not reach the global optimum [ 11 - 12 ], and there are 
three electable solutions exist. Generally, solutions 
with slightly better fitness (point C or B) prevailed, so 
the solution trapped into the valley prematurely. The 
more attractive solution(point G) is relatively far 
away from point A, but it nears the global optimal. To 
prevent prematurity, point G with slightly worse 
fitness than C, it needs a higher rank to be selected. 
That is, a higher RD k is awarded toa longer D k . 
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Stopping Criteria: Stopping criteria is given in the 
following order: 

1. Maximum allowable number of iterations reached. 

2. Number of iterations reached without improving 
the current best solution. 

Fig. 4 shows the flow chart of the proposed algorithm. 

C. Parameter identifications 

How to define the fitness function is the key point of 
the MPSO, since the fitness function is a figure of 
merit, and could be computed by using any domain 
knowledge. In this paper, we adopt the fitness 
function as follows [13-14]: 

n 

F {parameters) = y E? (28) 

7=1 

E. =x (0 -x* (0 (29) 

p 

where n is the total number of samples and j is the 

• *(z) 

calculated error of the 1 th sampling time, x is a 
solution by using the fourth-order Runge-Kutta 
method to solve the dynamic Equations (31) for the 
PM synchronous servomotor drive coupled with a 
slider-crank mechanism with the parameters identified 

(0 

from these two methods, and x is the displacement 
measured experimentally at the 1 th sampling time. 

IV. SIMULATIONS AND EXPERIMENTAL 
RESULTS 

A block diagram of the computer control system for 
the PM synchronous servomotor drive coupled with a 
slider-crank mechanism is shown in Fig. 5(a) and the 
experimental equipments of a slider-crank mechanism 
of computer control system is shown in Fig. 5(b). The 
control algorithm is implemented using a Pentium 
computer and the control software is LAB VIEW. The 
PM synchronous servomotor is implemented by 
MITSUBISHI HC-KFS43 series. The specifications 
are shown as follows: Rated output 400 (W), Rated 
torque 1.3 (Nm), Rated rotation speed 3000 (rpm) and 
Rated current 2.3 (A). The servo is implemented by 
MITSUBISHI MR-J2S-40A1. The control system is 
Sine-wave PWM control, which is a current control 
system. In order to measure the angle and angular 
speed of the disk and the position and velocity of the 
slider B, the interface of the device is implemented by 
motion control card PCI-7342. It can measure the 
angle of the disk and the position of slider B at the 
same time. 


The main parameters of a slider-crank mechanism and 
servomotor used in the numerical simulations and the 
experiments are as follows 
mj =0.232 kg, m 2 =0.332 kg, 
m 3 = 0.600 kg, r = 0.030 m, 
l = 0.217 m, F b = 0.100 TV, 

F e = 0.000 N, i q =0.400A, 

K, =0.5652 Nm/A, J m = 6.700-10” 5 Nms 2 , 

B m = 1.430-10“ 2 Nms!rad. 

The physical model of the slider-crank mechanism 
driven by a servomotor. In the parameter 
identification, we utilize the MPSO and PSO methods 
to identify the 5 parameters m 2 , m 3 , rand/ 

simultaneously, and the fitness function is described 
as Eq. (20). The identified results are given in Table 
1. 

Figure 6 shows the displacement errors of the slider- 
crank mechanism by both the PSO and MPSO 
methods. Furthermore, their displacement errors are 
all about ± 0.2 . Figure 7 shows the convergence 

characteristics in PSO and MPSO methods of the 
slider-crank mechanism system. It is seen from Fig. 7 
that the proposed MPSO method is superior to the 
PSO method. 

V. CONCLUSION 

The dynamic formulations of a slider-crank 
mechanism driven by a field-oriented PM 
synchronous motor have been successfully formulated 
with only one independent variable. Furthermore, the 
main objective of this study is to utilize PSO and 
MPSO methods to identify a slider-crank mechanism 
driven by a servomotor. According to the comparisons 
between identified results and displacement errors, it 
is found that MPSO method has the best matching 
with the experimental results. 

It is concluded that the implementations of MPSO are 
different from the PSO in five aspects. Firstly, its 
fitness function considers the distance to avoid 
converging to a local optimum. Secondly, for the 
MPSO, vectors with good enough fitness scores 
would be used as candidates to create new solutions. 
Thirdly, it has the advantage of the MPSO to conquer 
various constraints without using the fitness function 
with penalties, and can perform better. Fourthly, the 
solution is coded with a decimal representation, and 
saves computer memory. Finely, the gradually 
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decaying parameters can satisfy a successive statistic 
searching process by first using the diversification 
(bigger parameters) to reserve the larger attractive 
region. Then, the intensification (smaller parameters) 
used to search the small neighborhood of an elite 
solution. 
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Fig. 1(a) The physical model of a slider-crank 
mechanism. 



Fig. 1(b) The experimental equipments of a slider- 
crank mechanism 
Fig. 1 The slider-crank mechanism 


PM Synchronous Motor 
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Fig.2 Simplified control block diagram 



Fig.3 The concept of distances. 
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Fig.4 The flow chart of the MPSO. 



(a) 



(b) 

Fig.5 The experimental setup, (a) Computer control 
system block diagram; (b) The experimental 
equipments of a slider-crank mechanism of 


Computer control system. 

xlO -9 



Fig. 6 The comparisons in displacement errors 
between the PSO and MPSO methods. 



Fig. 7 Comparison of convergence characteristics 
in PSO and MPSO methods of the slider-crank 
mechanism system. 
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Table 1 

. The identifiec 

parameters of the numerical simulations. 

Parameter 

m x (kg) 

m 2 M 

(£g) 

r (m) 

l(m) 

Feasible domain 

0 . 000 - 1.000 

0 . 000 - 1.000 

0 . 000 - 1.000 

0 .000-0.100 

0 . 000 - 1.000 

The actual value 

0.232 

0.332 

0.600 

0.030 

0.217 

The identified value of the 
PSO method 

0.311 

0.304 

0.752 

0.024 

0.325 

The identified value of the 
MPSO method 

0.234 

0.331 

0.603 

0.030 

0.216 
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